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Abstract 

The often used GeU-Mann-Oakes-Renner (GMOR) mass formula for 
Nambu-Goldstone bosons in QCD, such as the pions, involves the con- 
densate <qq >, /tt and the quark current masses. Within the context 
of the Global Colour Model (GCM) for QCD a manifestly different for- 
mula was recently found by Cahill and Gunner. Remarkably Langfeld 
and Kettner have shown the two formulae to be equivalent. Here we ex- 
plain that the above recent analyses refer to the GCM constituent pion 
and not the exact GCM pion. Further, we suggest that the GMOR for- 
mula is generic. We generalise the Langfeld-Kettner identity to include 
the full response of the constituent quark propagators to the presence 
of a non-zero (and running) quark current mass. 
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1. Introduction 

The properties of the pion continue to be the subject of considerable theoretical 
and experimental interest in QCD studies. The pion is an (almost) massless Nambu- 
Goldstone (NG) boson and its properties are directly associated with dynamical chiral 
symmetry breaking and the underlying quark-gluon dynamics. Recently there has been 
renewed interest in the mass formula for the pion [Q, ^, Q and the relationship with the 
well known Gell-Mann-Oakes-Renner (GMOR) mass formula, as in (1) and (2). Here 
we extend the study of these relationships and show how one must carefully appreciate 
the different quantum field theoretic approaches that are actually being employed, often 
without explicit exposition. 

One expects that there should be some perturbative expression for the almost NG 

boson pion mass in terms of the small quark current masses which is built upon the 

underlying non-perturbative chiral-limit quark-gluon dynamics. While the relation of 

the low pion mass to the breaking of chiral symmetry dates back to the current algebra 

era and PCAC the often used implementation in QCD has the form, 

2 {mu + md)p 
= 72 

Jn 

where the integral p =<qq > is the so-called condensate parameter. For Nc = 3 

p = NMG{x = 0)) = 12 1 ^^siq'), (2) 

and /tt is the usual pion decay constant. (Note: our definition for p has an unconven- 
tional sign). In (2) (Ts{s) = crs(s; 0) is the chiral limit scalar part of the quark propagator 
which, utilizing only the Lorentz structure (and for later reference 0"„(s) = cj„(s;0)), 
we can write in full generality as 

G{q; m) = {iA{s] m)q.^ + B{s; m) + m{s))~^ = —iq.^ajj{s; m) + as{s; m), (3) 

from which we easily deduce that 

/ N ^ B{s;m) +m{s) 



A{s; m) 

sA{s] niY + {B{s; m) + m(s)y 



Here m{s) is the running quark current mass, but only the combination B{s; m) + m(s) 
appears. We note that the expression for p in (2) is divergent in QCD, because for 



s ^ oo B{s) decreases like l/{sln[s/A'^]^^^) where A = 12/(33 - 2Nf) and A is the 

QCD scale parameter. Some integration cutoff is usually introduced, and the values of 

m and <qq> are quoted as being relative to this cutoff momentum, often IGeV. The 

GMOR relation has been considered in various approaches, such as operator product 

expansions (OPE) HI, QCD sum rules |^, and recently, finite energy sum rules and 

Laplace sum rules . 

In a new mass formula for the pion mass was derived. The analysis in |Q] 

exploited the intricate interplay between the constituent pion Bethe-Salpeter equation 

(BSE) and the non-linear Dyson- Schwinger equation (DSE) for the constituent quarks, 

resulting in the new expression 
24771 r df'q 

"^^ = 'jr J J^:;^ ies{s)(Js{s) + s€y{s)a^{s)) c{s) + O(m^), (6) 
where c(s) is a naturally arising cutoff function 

= (7) 

Here es{s) and €y{s) are functions which specify the response of the constituent quark 
propagator to the turning on of the quark current mass; see ( p^ ) and (p3|). Note that 
the GMOR mass formula (||) and (^) appear to be manifestly different from the new 
expression in (^). However Langfeld and Kettner [^] have shown, by further analysis 
of the DSE for the constituent quark propagator, and ignoring for simplicity the e^(s) 
vector response term, that the two mass formulae are equivalent, even though the 
integrands are indeed different. 

Here we first demonstrate that in the quantum field theoretic analyses different 
concepts are often being used and confused in the literature. In this respect we carefully 
distinguish between the constituent pion and the full or exact pion, and their relevant 
mass expressions. Little detailed progress has been made in the exact analysis of QCD, 
and so we use the Global Colour Model (GCM) to illustrate these differences. Further 
we extend the Langfeld-Kettner identity [Q] to include the vector response function 
et,(s) and a quark runnning current mass function m(s). We show that the new mass 
formula can indeed be written in the GMOR form, with both (||) and (||) each now 
generalised to include a running current mass. 

To be clear we note that this report contains no analysis of the mass formula for 
the full pion in QCD, or even in the GCM. However if the GMOR relation is also the 



correct QCD result up to 0{m), then that would indicate that the GMOR relation is 
in fact a generic form that arises whether we are dealing with the full pion or with the 
constiutent pion, and whether we are analysing QCD or some approximation scheme 
to QCD, such as the GCM, provided we carefully preserve the dynamical consequences 
of the dynamical breaking of chiral symmetry and its activation by the underlying 
quark-gluon dynamics. We also note that the GMOR formlua has been 'derived' in the 
past, but such analyses in general appear to have brushed over the various subtleties 
presented herein. Ref.|§], in its appendix, illustrated this by using one example of an 
incorrect derivation leading to the GMOR relation. 

We exploit the GCM of QCD which has proven to be remarkably successful in 
modelling low energy QCD, as discussed in sec. 2. It should be noted that the GCM 
generates constituent hadrons which necessarily have the form of ladder states. The 
non-ladder diagrams then arise from hadronic functional integrations over constituent 
ladder hadrons. In sec. 3 we show the difference between the full or exact and the 
constiutent pion. This involves the use of effective actions and the fact that these 
effective actions refer to constituent hadrons. In sec. 4 the effective action for the chiral 
limit constituent pion is discussed. In sec. 5 the constituent quark propagators are 
given. They arise as the Euler Lagrange equations of the hadronised effective action 
for the GCM. They define the constituent quarks. Fluctuations about the minimum 
action configuration introduce constituent mesons, and these are described by ladder 
BSEs. Ad hoc alterations to these equations can introduce double counting problems. 
The full (observable) states are produced by dressing each of the constituent states by 
other states, as is made clear by the functional integral formalism in sec. 2. In sec. 6 
the constituent pion mass formula, (|^, is derived, but here generalised to include the 
quark running mass. In sec. 7 the Langfeld-Kettner identity is generalised to include 
the vector reponse function and a quark running current mass. This identity leads from 
the new formula in to the GMOR formula in (|^). 

2. The Global Colour Model of QCD 

An overview and an insight into the nature of the non-perturbative low energy 
hadronic regime of QCD is provided by the functional integral hadronisation of QCD 
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Q. In the functional integral approach correlation functions of QCD are defined by 

g{..,x,...)= j VqDqVA....q{x) exp{-SQCD[A,q,q]) (8) 

the kernel of which includes (not shown) gluon string structures that render Q{..,x, ...) 
gauge invariant. One example of (|8|) would be the pion correlation function, which may 
be defined by the connected part of 

Q^{x,y,z,w) = j VqVqVA....q{x)i^^nq{y)..q{z)t^^Tiq{w). .exp{-SQCD[A,q,q]). (9) 

The vrvr scattering amplitude, for example, is also defined by such a functional inte- 
gral. The pion mass is defined by the position of the pole, wrt the centre-of-mass 
(cm) momentum, of the Fourier transform of the translation invariant amplitude Q-,^. 
Apart from lattice computations a direct computation of these functional integrals is 
not attempted. Amplitudes, such as (^, when the on-mass-shell conditions are im- 
posed, define the observables of QCD, such as the pion. Theoretical analysis of these 
amplitudes proceeds by more circumspect techniques, some of which we clarify here. 

The correlation functions, as in @, may be extracted from the generating func- 
tional of QCD, ZQC'£)[7y, ry, ..], defined in (|lO|). However the interactions of low energy 
hadronic physics, such as vrvr scattering, are known to be well described by effective 
actions which refer only to hadronic states, although the various parameters in these 
effective actions could only be determined by fitting to experimental data. Hence we 
expect that the functional integrals, such as (|^, should also be extractable from a 
hadronic functional integral, as is indeed possible 

Zqcd^.V^A = j VqDqVAexp{-SQCD[A,q,q] +fiq + qr]) (10) 
« lv7rVNVN...exp{-Shad[7r,.;N,N,..] + J^[r],rj]TT + ..) (11) 

= Zhad[JAv,r]],-], (12) 

which produces a hadronic generating functional, Zhad[Jn\in,v]j --^^ ™ which source 
terms for the various hadrons are naturally induced. A partial derivation of this 
functional transformation proceeds as follows. First, and not showing source terms for 
convenience, the gluon integrations are formally performed (ghosts also not shown) 

VqDqVAexp{-SQCD [A, q, q] ) 



= J VqVqexp{— j q{'^.d + M)q+ 

where = g(x)^7^(7(x), and D^y{x) is the exact pure gluon propagator 

D^,{x-y)= j VAAl{x)Al{y)exp{-SQCD[A,m)- (14) 



A variety of techniques for computing D^y{x) exist, such as the ghionic DSE ||10|] and 
lattice simulations The terms of higher order than quartic in the quark fields 

are beyond our ability to retain in the analysis. Dropping these terms, beginning 
with D'^'^p, .., defines the GCM. Apart from lattice modellings of QCD no theoretical 
analyses have incorporated these higher order correlation functions in hadron studies. 
The remarkable success of the GCM suggests, for unknown reasons, that these terms 
are ineffective, at least in (colour singlet) hadronic states. Of course -D^^p, .• play an 
important role in the DSE approach [|lO| to estimating (]l^). This GCM truncation is 
equivalent to using a quark-gluon field theory with the action 

SGCM[q.q,A;] = J (^q(j.d + M{id)+iA';,^j,)q+U';,D-,Hid)A'^^ . (15) 

Here D^^ijp) is the matrix inverse of D^^{p), which in turn is the Fourier transform of 
Dfjiy{x). This action has only a global colour symmetry, unlike the local colour symme- 
try that characterises QCD. The gluon self- interactions, which arise as a consequence of 
the local colour symmetry in (14), lead to D~l{p) being non-quadratic. Its precise form 



is unknown, but see @ and 0. In the GCM a general non-quadratic form for D^i, is 
retained, modelling this significant property of QCD. One can even attempt to extract 
Dfiuip) from meson data [|l^]. The quark current masses M{s) = {mu{s) , md{s) , ..} in 
(^) are allowed to run. In momentum space this leads to the s-dependence of m{s) in 

Hadronisation involves a sequence of functional integral calculus changes of 
variables involving, in part, the transformation to bilocal meson and diquark fields, 
and then to the usual local meson and baryon fields (sources not shown). 



Z = J VqVqVAexpi-SQCDiA, q, q]+m + QV) 
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. / VirD,VAe.p{-SacM [A, g, g]+m + m) (GCM truncation) 

= J DBDVDV''exp{-Sbi[I3,V,V*]) (bilocal fields) (16) 

= JvTT...VNVN...expi-Shad[7^,-,N,N,..]) (local fields) . (17) 

The derived hadronic action that finally emerges from this action sequencing, to low 
order in fields and derivatives, has the form 

Shad[TT,...,N,N,..] = 

j d'^xtriNi-f.d + mN + ArriN - m7vV2i75vr"T" + ..)N} + 



+ / d^x 



f2 f2 



2 LV p . . TT J . 2 

J L 

2 

ifujfpf-KGi^pT^(i^i/(jr^pdvPu'dT'^~\~ 



+ -^e^^^rtriTT.Fd^TT.Fd^TT.Fd^TT.FdrTT.F) + 



(18) 



The bilocal fields in (|1^) arise naturally and correspond to the fact that, for in- 
stance, mesons are extended states. In (|9|) we can see that the pion arises as a correlation 
function for two bilinear quark structures. This bosonisation/hadronisation arises by 
functional integral calculus changes of variables that are induced by generalised Fierz 
transformations that emerge from the colour, spin and flavour structure of QCD. 

The final functional integration in ( |T7|) over the hadrons give the hadronic observ- 
ables, and amounts to dressing each hadron by, mainly, lighter mesons. The basic 
insight is that the quark-gluon dynamics, in (Ph, is fluctuation dominated, whereas 



the hadronic functional integrations in (17) are not, and for example the meson dress- 



ing of bare hadrons is known to be almost perturbative. In performing the change 
of variables essentially normal mode techniques are used Q. In practice this requires 
detailed numerical computation of the gluon propagator, quark propagators, and me- 
son and baryon propagators. The mass-shell states of the latter are determined by 
covariant Bethe-Salpeter and Faddeev equations. The Faddeev computations are made 
feasible by using the diquark correlation propagators; the diquarks being quark-quark 
correlations within baryons. 
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3. Constituent Hadrons 

We now come to one of the main points. Using the functional hadronisation we 
can write Qt^ in the form (0) or, from ([l7|), in the form: 



g^{X, Y)= J VTr..VNVN...7r{X)TT{Y)exp{-Shad[7T, A^, N, ..]), (19) 
in which X = and Y = are cm coordinates for the pion. We note that now 



the pion field appears in Sfiad[T^, X , N, ..] which is in the exponent of (19), and it 
appears with an effective-action mass parameter m-j^. As we now discuss, it is important 
to clearly distinguish between this mass (and the equations which define its value) and 
the pion mass that would emerge from (^ or (|l9|). Eqns.(|9|) or (|l9| ) define the observable 
pion mass. Whereas the mass in ( |T8|) defines the constituent pion mass. There is no 
reason for these to be equal in magnitude, though they may well both be given by the 
generic GMOR form. 

How do the constituent hadrons arise in (0)? In going from (|l^) to ([T7| ) an 
expansion about the minimum of St,i[B,'D,'D*] is performed. First the minimum is 
defined by Euler-Lagrange equations (ELE) 

— = and — = 0. (20) 



These equations have solutions B ^ and D = 0. Eqns.(20) with P = is seen, after 
some analysis, to be nothing more than the DSE for the constituent quark propagator in 
the rainbow approximation (see (|27| ) and (^) in sec. 5). The occurrence of the rainbow 
form of these equations is not an approximation within the GCM. The non-rainbow 
diagrams, corresponding to various more complicated gluon dressing of the quarks, are 
generated by the additional functional integrals in ([T7|). Ad hoc alterations to these DSE 
constituent quark equations will lead to double counting of certain classes of diagrams. 
The generation of a minimum with ;B 7^ is called the formation of a condensate, here 
a qq condensate. That D = means that in the GCM no diquark or anti-diquark type 
condensates are formed. 

Next in going from ( [l6| ) to ([T7| ) we must consider the fluctuations or curvatures 
of the action for the bifocal fields. One finds that the curvature 6'^Shi/SB6B when 
inverted gives the meson propagators, but with only ladder gluon exchanges. Again 
non-ladder diagrams are generated by the functional hadronic integrations in ([T7|). 



Similarly inversion of the curvatures in the diquark sector 5'^ Si,i/5'D6'D* leads to diquark 
propagators, but with only ladder gluon exchanges between the constituent quarks. 

We note that the generalised bosonisation with meson and diquark fields leads 
to some additional complications that we shall consider elsewhere, but which do not 
impinge on the basic point being made here. This meson-baryon hadronisation is based 
upon a generalised Fierz transformation that induces the appropriate colour singlet 
anti-quark - quark correlations, and colour anti-triplet quark-quark correlations that are 
in the correct colour state for quark-quark correlations within a colour singlet baryon. 



An earlier bosonisation [13| used a Fierz transformation that lead to only the 
meson sector of the GCM. In this bosonisation the meson effective action involves con- 
stituent states that are generated by naturally arising and exclusively rainbow or ladder 
diagrams. Then all of the other diagrams contributing to the observable states are gen- 
erated by the functional integrations in (|T7|). Hence we see that in the exponent in 
(0) there arise particular propagators for quarks, mesons, diquarks and even baryons. 
These propagators and their associated fields will be defined to be the constituent 
states. They could also be described as core states. The observables are generated by 
the hadronic functional integrals in ([TtI) , and correspond to dressing each constituent or 
core state with other such states. Hence the hadronic effective action in (|l^ ) contains 
a variety of parameters that refer to the constituent states. 

Nevertheless one often compares these parameters with the parameter values for 
the fully dressed constituent states, that is the observable hadronic states. This ap- 
pears to be valid because in general the dressing produces only a small shift in the 
parameter values. However one known exception is the nucleon where pion dressing of 
the constituent nucleon state reduces its mass by some 200 — SOOMeV. This mass shift 
emerges from consideration of the functional integral 

gNiX,Y) = J VTT..NiX)N{Y)expi-Shad[TT,...,N,N,..]), (21) 

where mainly the pions, but as well other mesons are used to dress the nucleon. Of 
course one usually casts this into the form of a non-linear integral equation for the 
meson dressed nucleon correlation function. 

4. Chiral Limit Constituent Pions 
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When the quark current masses — > SqcdIi-, has an additional global 

Ul{Nf) ® Ur{Nf) chiral symmetry: writing q-f^q = qRjf.qR + ^lTm^l where g^.L = 
Pr,lQ and q^ l — qPl.r we see that these two parts are separately invariant under 
QR URqR,qji IrUr and qi Uiqi^QL ^ '^lU'I- Its consequences may be 
explicitly traced through the GCM hadronisation. First the ELEs 5Si,i/5B=0 have 
degenerate solutions. In terms of the constituent quark propagator this degeneracy 
manifests itself in the form 

Giq; V) = [iAiq)q.^ + VB{q)]-^ = C^G{q; 1)C^ (22) 
where 

C = VV, V = exp{i^2j5Tr^F^) (23) 

in which the {vr"} are arbitrary real constants. The degeneracy of the minimum implies 
that some fluctuations in 5'^ Sbi/ 5B5B have zero mass; these are the NG BSE states, 
and this indicates the realisation of the Goldstone theorem. 

In the hadronisation, in going from (ID to (0), new variables are forced upon 
us to describe the degenerate minima (vacuum manifold). This is accomplished by a 
coordinatisation of the angle variables {vr"}: 

U{x) =exp{i^2TT''{x)F'') (24) 

V{x) = PlU{x)^ + PrU{x) = exp(V275^"(a;)F") (25) 

The NG part of the hadronisation then gives rise to the constituent pion effective 
action 

j d\ (^ltr{d^Ud^U^) + KMd^Ud^U^) + |tr([l - ^^±^]7W)+ 

+K2tr{[d^Ud^U^]^) + K^trid^UdM^d^UdM^) + (26) 



This is the ChPT effective action |14|j , but with the added insight that all coefficients 
are given by explicit and convergent integrals in terms of A and B, which are in turn 
determined by D^^. The higher order terms contribute to vrvr scattering. The depen- 



dence of the ChPT coefficients upon D/^i, has been studied in [|l5| , 12, 16 1 in which the 



GCM constituent pion expressions for the various parameters were used. However in 
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view of the apparent generic role of the GMOR relation one should keep in mind the 
possibility that the functional form of the dependences of the parameters ki, .. on the 
quark correlation functions A and B might also be generic. At present the final func- 
tional integral dressing to obtain the pion observables has not been carried out. This 
amounts to the assumption that the constituent pion forms are sufficiently accurate. 
The hadronisation procedure also gives a full account of NG-meson - nucleon coupling. 

The GCM is in turn easily related to a number of the more phenomenological 
models of QCD, as indicated in fig.l. They include the Nambu-Jona-Lasinio Model 
(NJL) 0, ChPT ig, MIT and Cloudy Bag Model (CBM) [|l§, Soliton Models 0, 
Quantum Hadrodynamics (QHD) and Quantum Meson Coupling model (QMC) 
II21II . We also indicate that the pure gluon correlation function in (|l^ ) may be obtained 
from lattice computations and used in the GCM. The relationships indicated in fig.l 



are discussed in [22|. 



QCD 



GCM 




Lattice Gluons 



NJL 



ChPT 



MIT, Cloudy Bag, 

Soliton Models, QHD, QMC. 

Lattice Hadrons 



Fig.l. Relationship of the GCM to QCD and other models 



5. Action Minimum and Pionic Fluctuations 

The GCM involves the solution of various integral equations for the constituent 
correlation functions. As we saw in sec. 3, the first equation involves the determination 
of the minimum of the bilocal effective action, and this reduces to solving the DSE for 
the constituent quark propagator necessarily in the rainbow form (the so-called vacuum 
equation of the GCM g 

Rr 2 , 16 /■ A p. ^ B{q^;m)+m{q'') 

^[P ,rn) 3 y (27r)4'^^^ ^ ' q^ Aiq^ ; mf + {B {q^ ; m) + m{q^))^ ' ^ ' 

T A, 2 N -.12 8 /■ (fq ^, , A{q^]m) 

[A{p ;m -Ip = - / ——^q.pD{p-q). — r^rrf, 28) 
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For simplicity we use a Feynman-like gauge in which Di^i^{p) = 5^uD{p) (the quark- 
gluon couphng is incorporated into D). The formal results of the analysis here are not 
gauge dependent. Even in numerical studies the Landau gauge can also be used; see 
p^ . We have also included, for generality, a running current mass for the quarks. 
Using Fourier transforms ([27| ) may be written in the form, here for m = 0, 

16 as{x) 

which implies that knowledge of the quark propagator determines the effective GCM 
gluon propagator. Multiplying (p9| ) by B{x)/ D{x), and using Parseval's identity for 
the RHS, we obtain the identity 

d^x-^ = - / —^,B{q)as{q). (30) 



D{x) 3 J (27r)4^ 
The second basic equation is the ladder form BSE for the constituent pion mass- 
shell state, which arises from the mesonic fluctuations about the minimum determined 



by ( [27| ) and (28). Again this ladder form cannot be generalised without causing double 
counting of some classes of diagrams at a later stage, and without also damaging the 
intricate interplay between (p7|) , (p8|) and the BSE 

Tf{p,P) = ^D{p-q)trsFiG+T^G^Tf)THq,P) (31) 

where G± = G{q it y). This BSE is for isovector NG bosons, and only the dominant 
r = T^THj^ amplitude is retained (see |23] for discussion); the spin trace arises from 



projecting onto this dominant amplitude. Here {T^, b = 1, .., Np — 1} are the generators 
of SU{Nf), with tr{TfT9) = \6fg. 

The BSE ( |3l|) is an implicit equation for the mass shell = — M^. It has solutions 
only in the time-like region P^ < 0. Fundamentally this is ensured by (^^ and (|2^ ) 
being the specification of an absolute minimum of an effective action after a bosonisation 
Q. Nevertheless the loop momentum is kept in the space- like region q^ > 0; this mixed 
metric device ensures that the quark and gluon propagators remain close to the real 
space-like region where they have been most thoroughly studied. Very little is known 
about these propagators in the time-like region q^ < 0. 

The non-perturbative quark-gluon dynamics is expressed here in (|27| ) and (|2^ . 
Even when m = ( ^7|) can have non-perturbative solutions with B ^ 0. This is the 
dynamical breaking of chiral symmetry. 
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When m = (^) has a solution for = 0; the Goldstone theorem effect. For 
the zero Hnear momentum state {Pq = 0, P = 0} it is easily seen that (^) reduces to 
{^) with Tf{q,0) = B{q^). When P / then Tf{q,P) / B{q), and (|3l|) must be 
solved for Tf{q,P). 

6. Constituent Pion Mass 

We shall now determine an accurate expression for the mass of the constituent pion 
when m{s) is small but non-zero. This amounts to finding an analytic solution to the 



BSE (13 1|) , when the constituent quark propagators are determined by ( \27\ j and (|28D . 
The result will be accurate to order m. 

For small m 7^ we can introduce the Taylor expansions in m{s) 

B{s; m) + m(s) = B{s) + m{s).es{s) + 0{m'^), (32) 

A{s; m) = A{s) + m(s).e„(s) + ©(m^). (33) 

For large space-like s we find that ^ 1, but for small s we find that es{s) can be 
significantly larger than 1. This is an infrared region dynamical enhancement of the 
quark current mass by gluon dressing, and indicates the strong response of the chiral 
limit constituent quark propagator to the turning on of the current mass. A plot of es{s) 



is shown in Q. Higashijima |24| and Elias [25| have also reported similar enhancements 



of the current quark masses in the infrared region. 

Even in the chiral limit the constituent quark running mass M(s) = B{s)/A{s) is 
essential for understanding any non-perturbative QCD quark effects. The integrand of 
a BSE contains the gluon correlation function, constituent quark correlation functions 
and the form factor for the state (see for example (pl|)). This integrand shows strong 



peaking at typically s « 0.3^6^^. At this value we find |T2| that M{s) ^ 270MeV. 
This is a property of the constituent hadrons. It does not include any effects from the 
dressing of these hadrons via (^). This mass is called the constituent quark mass. 
Because of the infrared region enhancement of the quark current mass we find that this 
constituent mass rises quickly with quark current mass; see |12| . 

Because the pion mass tjItj is small when m is small, we can perform an expansion 
of the dependence in the kernel of (|3l|). Since the analysis is Lorentz covariant we 
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can, without loss of validity, choose to work in the rest frame with P = {iraj^, 0) giving, 
for equal mass quarks for simplicity 

2 ^ r 



Tip) = -ml J ^D{p-q)I{snq) + 



(2vr) 

+y / j2^^^^~'^K{A{s) + m{s).e,is)y + {B{s) + m{s).esis)y^^'^^^---^^ ^^^^ 

where 

I{s) = 6 [a^ - 2{a,a', + sa^a',) - sia^a'^ - (a^f) - s^a^a'^ - (a'J)) . (35) 
By using Fourier transforms the integral equation (jsj), now with explicit depen- 



dence on m^r, can be expressed in the form of a variational mass functional, 

^ U[T]^J i27Tys{A{s) + m{s).e,is)r + {Bis) + m{s).es{sW 

+ ^—fd'x^ (36) 

in which 

Urf = I ^I{s)T{q)'. (37) 



. (27r)4 

The functional derivative 5m^[r]^/(5r(g) = reproduces (|34[). The mass functional 
(^) and its minimisation is equivalent to the constituent pion BSE in the near chiral 
limit. To find an estimate for the minimum we need only note that the change in 
from its chiral limit value of zero will be of 1st order in m, while the change in the zero 
linear momentum frame T{q) from its chiral limit value B{q^) will be of 2nd order in 
m. 

Hence to obtain to lowest order in m, we may replace T{q) by B[q^) in 
and we have that the constituent pion mass is given by 
2 24 [ d^q Ms)B{s) + se,{s)A{s) B{s)^ 



U[BY J (27r)4 s^(s)2 + S(s)2 syl(s)2 + 5(s)2 

^ ""^'^^ ^^fd^^^^Oim^) (38) 



U[B]^J {27t)^ sAis)^ + B{s)^ 2U[BYJ D{x) 

However the pion mass has been shown to be zero in the chiral limit. This is confirmed 
as the two 0{mP) terms in (|3^ ) cancel because of the identity (|30|). Note that it 
might appear that fn would contribute an extra m dependence from its kernel in 
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However because the numerator in ( |3q ) is already of order m, this extra contribution 
must be of higher order in m. 

Hence we finahy arrive at the analytic expression, to 0{m), for the constituent NG 
boson (mass)^ from the solution of the BSE in (^) 

9 24 f d^q , fAs)B{s)^seAs)A(s) Bisf 9^ 



Eqn.(^) or ( p9D is the new form of the NG mass formula derived in [Q]. It would 
appear that expression (|6|) is manifestly different to the conventional GMOR form in 
(H) and @. However in the next section we generalise an identity found by Langfeld 
and Kettner Q which shows these forms to be equivalent. 

7. Relating the Mass Formulae 

Inserting (^) and (|3^ ) into (p7|), and expanding in powers of m(s), we obtain up 
to terms linear in m, and after using (p7|) with m = to eliminate the 0{mP) terms, 



m{p )es{p )=m{p) + -j ^^D{p - ^)^J^^yr^BW? 



3 y (27r)4 ^'(g2A(g2)2 + 5(^2)2)2 

1« f "'^Dtp-ri ^'^^'-^'^'^-'g^^y) . (40) 



3 7 (27r)4 (g2A(g2)2 + 5(^2)2)2 

We now multiply throughout by i?(p^)/(p^74(p^)^ + i?(p^)^), and integrate wrt 
p. Using again the chiral limit of (27) there is some cancellation of terms, and we are 



left with a generalised Langfeld-Kettner identity 



-S(p^)^ / B{p^)m{p^)es{p'^) p^ A{p^)m{p^)e^{p^) 



p^A{p^)'^ + B(p2)2 \^p2^|-^2)2 _^ _B(p2)2 p2^(^p2~)2 _^ ^(^2)2 j 

, m{j?)B(j?) 

^ p^Aip^Y + B(p2)2 ^^^^ 

Remarkably, noting (Q) and (|5|), we see that using this identity in (^) or (|3^) finally 
completes the derivation of the GMOR expression for the mass of the constituent pion. 

We thus see that despite its apparently simple form the GMOR expression actually 
depends on two identities that follow from the non-linear constituent quark DSE, as 
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well as on the subtle interplay between this constituent quark equation and the BSE for 
the constituent pion. These in turn arise from the careful self-consistency rendered by 
the functional integral prescription which ensures that the fluctuation spectrum for the 
bilocal action is precisely related to the Euler-Lagrange equations. Ad hoc alterations 
will invalidate this connection and therefore the derivation of the GMOR expression. 

8. Conclusion 

We have indicated the careful considerations that must be given to modelling QCD 
via the GCM and the manner in which this leads to hadronic effective actions. We have 
defined constituent quarks, meson, diquarks and baryons as those states that appear 



in the effective action, i.e. in the exponent, as in (17). These constituent states are 
then further dressed by the functional integrations in ([T7|). Remarkably this GCM 
structuring of the quantum field theoretic analysis implies, at least in the simplest 
version of the GCM, that the constituent states are described by sums of rainbow or 
ladder diagrams, and that the functional integrations then build up all the remaining 
diagrams, amounting to the vast array of crossed diagrams and vertex dressings etc. 
Because in most cases these extra dressings do not cause large changes in the values of 
the constituent masses, coupling constants,., the effect of the inclusion of these extra 
diagrams is not manifestly large. Of course this is not surprising because the GCM 
hadronisation allows us to assess the significance of a constituent state through its mass; 
low mass states should be more important than very massive states. This implies that 
the pion dressing is the largest such effect. However inclusion of this dressing for the 
constituent nucleon is known to be significant, and is a result of the inclusion of various 
non-ladder diagrams in the observable nucleon. 

We have also carefully indicated that it is the mass of the constituent pion that is 
analysed here and, by using various identities that follow from the non-linear equation 
for the constituent quark, one can show that the mass of this constituent pion is indeed 
given by the GMOR formula, with the scalar part of the constituent quark correlation 
function appearing. This does not preclude the fact that presumably the observable 
pion also has its mass obeying a GMOR formula, but one in which the full quark scalar 
correlation function appears. That is, the GMOR relation is generic. 
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